We perform a time-dependent Dalitz plot analysis of Kgir^-K^ in 

order to extract the CP violation parameters of fo{980)Kg and p^{770)Kg 
and direct CP-asymmetries of i^*^(892)7r~. The results are obtained from a 
data sample of (383 it 3) x 10^ BB decays, collected with the BaBar detector 
at the PEP2 asymmetric-energy B factory at SLAC. The measured values 
of 2/3eff in fiO decays to /o(980)E:0 pO{770)K^ are (89^^^ ± 5 ± 8)° and 
(37]']}y ib5zb6)°, respectively, where the first quoted uncertainty is statistical, 
the second is systematic and the third is Dalitz plot signal model uncertainty. 
We measure the significance of 2/?efT(/o(980)i^g) 7^ to be 4.3a. In decays 
to i^*(892)7r we find Acp = -0.18 ± 0.10 ± 0.03 ± 0.03. The measured 
phase difference between the decay amplitudes of B^ i('*+(892)7r~ and 
S° K*-{892)tt+ is (-164 ± 24 ± 12 ± 15)°. All results are preliminary. 
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1 INTRODUCTION 



The Cabibbo-Kobayashi-Maskawa (CKM) mechanism [1, 2] for quark mixing describes all transi- 
tions between quarks in terms of only four parameters: three real rotation angles and one irreducible 
phase. Consequently, the flavor sector of the Standard Model (SM) is highly predictive. One partic- 
ularly interesting prediction is that mixing-induced CP asymmetries in decays governed hy h ^ qqs 
(g = n, s) transitions are, to a good approximation, the same as those found in 6 ^ ccs tran- 
sitions. Since flavor changing neutral currents are forbidden at tree-level in the Standard Model, 
the h s transition proceeds via loop diagrams (penguins) , which are affected by new particles in 
many extensions of the SM. 

Recently, various different b ^ s dominated charmless hadronic B decays have been studied in 
order to probe this prediction. The values of the mixing-induced CP asymmetry measured for each 
(quasi-)two-body mode can be compared to that measured in 6 — > ccs transitions (typically using 

— J/il:Kg). A recent compilation [3] of these results shows that they tend to have central values 
below that for b — ^ ccs. However, there is currently no convincing evidence for new physics effects 
in these transitions. The most recent theoretical evaluations [4-12] suggest that SM corrections to 
the b — > qqs mixing-induced CP violation parameters should be small (in particular for the modes 
rj'K^ and K^K^Kg), and tend to increase the values (i.e. the opposite trend to that seen in 
the data). Clearly, more precise experimental results are required. 

The compilation given in [3] includes several three-body modes, which may be used either 
by virtue of being CP eigenstates [KgKgKg, Kg-K^n^) [13] or since their CP content can be 
determined experimentally {K'^K~ K^) [14, 15]. It also includes quasi-two-body (Q2B) modes, 
such as /o(980)i<'^ and p^{770)Kg, which arc reconstructed via three-body final states (i^^^7r+7r~ for 
these modes). For these channels, the precision of the Q2B approach is limited as other structures 
in the phase space may cause interference with the resonances considered as signal. Therefore, more 
precise results can be obtained using a full, time-dependent Dalitz plot (DP) fit of 5° — ^ Kgir'^Tr^. 
Furthermore the interference terms allow the cosine of the effective weak phase difference in mixing 
and decay to be determined, helping to resolve ambiguities which arise from the Q2B analysis. 
This approach has recently been successfully used in a time-dependent Dalitz plot analysis of 
5° ^ K+K-K^ [15]. 

The discussion above assumes that the b ^ s penguin amplitTidc dominates the decay. However, 
for each mode contributing to the K^-k^tt^ final state, there is also the possibility of a 6 ^ u tree 
diagram. These are doubly CKM suppressed compared to the b ^ s penguin (the tree is O(A^) 
whereas the penguin is C>(A^), where A is the usual Wolfenstein parameter [16, 17]). However, 
hadronic factors may result in a relative enhancement and hence significant "tree pollution" . The 
relative magnitudes of the tree and penguin amplitudes, \T/P\, thus can be different for each 
Q2B state, as can the strong phase difference. However, the relative weak phase between them 
is of course the same - and in the Standard Model is approximately equal to 7. An amplitude 
analysis, in contrast to a Q2B analysis, yields sufficient information to extract relative phases and 
magnitudes. DP analysis of B^ KgTT^TT^ (and similar modes) can therefore be used to determine 
7 [18-21]. A comparison of the value obtained with that extracted from tree-level B DK decays 
provides a SM test. 

No results on time-dependent DP analysis of B^ — if^Tr+Tr" have yet been published. Belle 
have presented results of an analysis using 388 million BB pairs [22], which does not take into 
account either time-dependence or flavour-tag dependence. The results of the Belle analysis are 
consistent with other studies of the B^ — > i^^Tr+Tr" decay [14, 23], as well as with measurements 
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obtained from other Kirir systems: K'^tt'tt'^ [24, 25] and K+tt+tt" [26, 27]. The latter results 
indicate evidence for direct CP violation in the p^K'^ channel. If confirmed, this will be 

the first observation of CP violation in the decay of any charged particle. Taken together with the 
observation of direct CP violation in —>■ K~^ir~ decays [28, 29], these results suggest that large 
CP violation effects are possible in — > K*~^Tr~ (although current measurements of the effect are 
consistent with zero [23]). 

In this paper we present preliminary results from the first time-dependent Dalitz plot analysis 
of the B^ Kg7r^7T~ decay. In Section 2 we describe the time-dependent Dalitz plot formalism, 
and introduce the signal parameters that are extracted in our fit to data. In Section 3 we briefly 
describe the BaBMI detector and the data set. In Section 4, we explain the selection requirements 
used to obtain our signal candidates and suppress backgrounds. In the same section we describe 
the methods used to control experimental effects such as resolution in the fit to data. In Section 5 
we present the results of the fit, and extract parameters relevant to the contributing Q2B decays. 
In Section 6 we discuss systematic uncertainties in our results, and finally we summarize our results 
in Section 7. 



2 ANALYSIS OVERVIEW 

Using a maximum- likelihood fit, we measure relative phases and magnitudes for the different res- 
onant decay modes, taking advantage of the interference between them in the Dalitz plot. Below 
we detail the formalism used in the present analysis. 

2.1 DECAY AMPLITUDES 

We consider the decay of a spin-zero B^ with four-momentum ps into the three daughters tt"*", tt" 
and Kg with p_ and po, their corresponding four- momenta. Using as independent (Mandel- 
stam) variables the invariant squared masses 

s+ = {p++Pof , s- = {p-+pof , (1) 

the invariant squared mass of the positive and negative pion, sq = {p+ + p-)"^, is obtained from 
energy and momentum conservation 



So = rn^BO + + m^K° - s+ - s- . (2) 



The differential B^ decay width with respect to the variables defined in Equation (1) (i.e. the Dalitz 
plot) reads 

dr(B»^A>+.-) = (3) 

where A is the Lorentz-invariant amplitude of the three-body decay. In the following, the ampli- 
tudes A and its CP conjugate A correspond to the transitions B^ KgTT^7r~ and B^ — > i^^7r'''7r~, 
respectively. We describe the distribution of signal events in the DP using an isobar approximation, 
which models the total amplitude as resulting from a sum of amplitudes from the individual decay 
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channels 

N 

As+,s-) = ^CjFj{s+,S-) (4) 

TV 

A{s+,s-) = ^CjFj{s+,S-) (5) 
i=i 

where Fj are DP dependent dynamical amplitudes described in the following, and Cj are complex 
coefHcients describing the relative magnitude and phase of the different decay channels. All the 
weak phase dependence is contained in Cj, and Fj contains strong dynamics only, therefore 

Fj{s+,s.) = Fj{s.,s+). (6) 

The resonance dynamics are contained within the Fj terms, which are represented by the product 
of the invariant mass and angular distribution probabilities, i.e. 

if s_) = Rj{m) X r) x r) x T,(L,p, q) , (7) 

where 

• m is the invariant mass of the decay products of the resonance, 

• Rj{m) is the resonance mass term or "lineshape" (e.g. Breit-Wigner), 

• Xl are Blatt-Weisskopf barrier factors [30] with parameter r. These factors are taken to be 
unity in the present analysis, both in Equation (7) and in the lineshapes. The effect of this 
choice is accounted for systematic uncertainty. 

• p* is the momentum of the bachelor particle evaluated in the rest frame of the B, 

• p and q arc the momenta of the bachelor particle and one of the resonance daughters respec- 
tively, both evaluated in the rest frame of the resonance, 

• L is the orbital angular momentum between the resonance and the bachelor, and 

• Tj{L,p,q) is the angular distribution, where: 

L = : Tj = l , (8) 
L = l : Tj = -Ap-q, (9) 

L = 2 : Tj = l[3{p.q)'-{\p\\q\?] ■ (10) 

The lineshape differs for each component included in the fit. The lineshapes used are Relativistic 
Breit-Wigner (RBW) [31], Flatte [32], Gounaris-Sakurai (GS) [33] and LASS [26, 34, 35]. A 
flat phase space term has been included in the signal model to account for nonresonant (NR) 

—>■ KgTr~^Tr~ decays. The components of the signal model are summarized in Table 1. 
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Resonance 


Parameters 


Form Factor 


Rcf. for 
Parameters 


/o 


mass = 965 it 10 
5^ = 165 ± 18 
gK = 695 ± 93 


Flatte 


[36] 


A 


mass = 775.5 ± 0.4 
width = 146.4 ± 1.1 


GS 


[31] 


K*+(892) 


mass = 891.66 ± 0.26 


RBW 


[31] 


K* (892) 


width = 50.8 ± 0.9 






iir*+ (1430) 


mass = 1415 ± 3 


LASS 


[26, 37] 


K*-(1430) 


width = 300 ± 6 
cutoff = 2000 
a = 2.07±0.1(GeV-^) 
r = 3.32 zb 0.34 (GeV ) 






/x(1300) 


mass = 1449 dr 13 
width = 126 ± 25 


RBW 


[22] 


fof 127(11 


mass — 1 975 4 + 1 1 
width = 185.2l|^ 


RRW 




Xco(lP) 


mass = 3414.75 ± 0.35 

width = 10.1 ± 0.7 


RBW 


[31] 


NR decays 




flat phase space 





Table 1: Parameters of the DP model used in the fit. Values are given in MeV, unless mentioned 
otherwise. 

2.2 TIME DEPENDENCE 

With A.t = isig — ^tag defined as the proper time interval between the decay of the fully reconstructed 
— K^tt'^'k' (i?sig) and that of the other meson [B^^^) from the T(45), the time-dependent decay 
rate |^^g(Ai)p {\A~:,JyM)\^) when the B^^^ is a 5° (5°) is given by 



g-|At|/TgO 



1^1' + 1^1' 



4rB0 

+ {\A\^ - cos(AmdAt) 
± 2Im [AA*] sin(AmdAt) 



(11) 



where r^o is the neutral B meson lifetime and Am,^ is the B^B^ mass difference. In the last formula 
and in the following, the DP dependence of amplitudes is implicit. Here, we have assumed that 
there is no CP violation in mixing, and have used a convention whereby the phase from B^B^ 
mixing is absorbed into the B^ decay amplitude ^ {i.e. into the Cj terms). Lifetime differences in 
the neutral B meson system are assumed to be negligible. 



^In other terms, we assume that \q/p\ = 1 and absorb qjp into Cj 
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2.3 THE SQUARE DALITZ PLOT 



Both the signal events and the combinatorial e+e~ QQ (q = u,d,s,c) continuum background 
events populate the kinematic boundaries of the DP due to the low final state masses compared 
with the mass. We find the representation in Equation (3) is inconvenient when one wants to 
use empirical reference shapes in a maximum-likelihood fit. Large variations occurring in small 
areas of the DP are very difficult to describe in detail. We therefore apply the transformation 

ds+ ds_ — I det J| dm' dO' , (12) 

which defines the Square Dalitz plot (SDP). The new coordinates are 

m' = - arccos ( 2 "^o " "^r" _ = 1^, , (13) 

where mo = -^/sq is the tt+tt" invariant mass, m™^ = m^o — m^o and m™"^ = 2m^+ are the 
kinematic limits of mo, ^o is the vr+vr" resonance helicity angle, defined as the angle between the 
7r~ and the in the ■k~^it~ rest frame, and J is the Jacobian of the transformation. Both variables 
range between and 1. The determinant of the Jacobian is given by 

|detJ| = 4|p+||po|mo--^-^^ , (14) 



where |p+| = E*^ - m^^ and |Pol = ^ - m^^o , and where the 7r+ {K^) energy E\ {Eq), is 

defined in the vr+vr" rest frame. This same transformation has been used in previous B decay DP 
analyses, e.g. Ref. 



3 THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the BABAR detector at the PEP-II asymmetric- 
energy e+e" storage ring at SLAC between October 1999 and August 2006. The sample consists of 
about 347 fb~\ corresponding to (383 ± 3) x 10^ BB pairs collected at the T(45) resonance ("on- 
resonance"), and an integrated luminosity of 36.6 fb~^ collected about 40 MeV below the T(4S') 
( "off-resonance" ) . 

A detailed description of the BABAR detector is presented in Ref. [39]. The tracking system 
used for track and vertex reconstruction has two components: a silicon vertex tracker (SVT) and a 
drift chamber (DCH), both operating within a 1.5 T magnetic field generated by a superconducting 
solenoidal magnet. Photons are identified in an electromagnetic calorimeter (EMC) surrounding 
a detector of internally reflected Cherenkov light (DIRC), which associates Cherenkov photons 
with tracks for particle identification (PID). Muon candidates are identified with the use of the 
instrumented flux return (IFR) of the solenoid. 



4 EVENT SELECTION AND BACKGROUND SUPPRESSION 

We reconstruct B^ — > i^^Tr+Tr" candidates from pairs of oppositely-charged tracks and a -^■ 
7r+7r~ candidate, which are required to form a good quality vertex. In order to ensure that all 
events are within the DP boundaries, we constrain the invariant mass of the flnal state to the B 
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mass. For the 7r"'"7r~ pair from the B, we use information from the tracking system, EMC, and 
DIRC to remove tracks consistent with electron, kaon and proton hypotheses. In addition, we 
require at least one track to be inconsistent with the muon hypothesis based on information from 
the IFR. The Kg candidate is required to have a mass within 15MeV/c^ of the nominal mass 
[31] and a decay vertex well separated from the decay vertex. In addition, the cosine of the 
angle between the Kg flight direction and the vector between the S-daughter pions and the Kg 
vertices must be greater than 0.999. 

A B-meson candidate is characterized kincmatically by the energy-substituted mass rriES = 



y (s/2 + Pi ■ pB)^/Ef — and energy difference AE = — ^\/s, where {Eb, Pb) and {Ei,pi) 
are the four-vectors of the i?-candidate and the initial electron-positron system, respectively. The 
asterisk denotes the T{4S) frame, and s is the square of the invariant mass of the electron-positron 
system. We require 5.272 < mES < 5.286 GeV/c^ and \AE\ < 0.065 GeV. 

Backgrounds arise primarily from random combinations in continuum events. To enhance dis- 
crimination between signal and continuum, we use a neural network (NN) [40] to combine four 
discriminating variables: the angles with respect to the beam axis of the B momentum and B 
thrust axis in the T{AS) frame; and the zeroth and second order monomials Lo,2 of the energy flow 
about the B thrust axis. The monomials are defined by L„ = pj x |cos 9i\^, where 9i is the angle 
with respect to the B thrust axis of track or neutral cluster i, pi is its momentum, and the sum 
excludes the B candidate. The NN is trained using off-resonance data as well as simulated signal 
events that passed the selection criteria. The final sample of signal candidates is selected with a 
requirement on the NN output that retains 90% (29%) of the signal (continuum). 

The time difference At is obtained from the measured distance between the z positions (along 
the beam direction) of the B^^^ and B^^^^ decay vertices, and the boost P'y = 0.56 of the e+e" 
system: At = Az/Pjc. B^ candidates with \At\ > 20 ps are rejected, as well as candidates for 
which the error on At is higher than 2.5 ps. To determine the flavor of the B^^^ we use the B 
flavor tagging algorithm of Ref. [41]. This produces six mutually exclusive tagging categories. We 
also retain untagged events in a seventh category to improve the efficiency of the signal selection 
and because these events contribute to the measurement of direct CP violation [42] . Events with 
multiple B candidates passing the full selection occur between ~ 1% of the time for nonresonant 
signal and ~ 8% of the time for B^ foidSO) Kg signal. If an event has more than one candidate, 
we select one using an arbitrary but reproducible procedure based on the event timestamp. 

With the above selection criteria, we obtain a signal efficiency determined from Monte Carlo 
(MC) simulation of 21 — 25% depending on the composition of the DP. 

Of the selected signal events, 8% of B^ p^K^, 6% of B^ /i:*(892)+7r- and 4% of B^ 
fo{98Q)K^ events are misreconstructed. Misreconstructed events occur when a track from the 
tagging B is assigned to the reconstructed signal candidate. This occurs most often for low- 
momentum tracks and hence the misrcconstiTicted events are concentrated in the corners of the 
DP. Since these are also the areas where the low-mass resonances overlap strongly with other 
resonances, it is important to model the misreconstruced events correctly. The details of the model 
for misreconstructed events over the DP is detailed in Section 4.2.1. 



We use MC simulated events to study the background from other B decays. More than fifty channels 
were considered in preliminary studies, of which twenty are included in the final likelihood model 
- decays with at least two events expected after selection. These exclusive B background modes 




4.1 BACKGROUND FROM OTHER B DECAYS 



13 



Mode 




Varied 


BR 


Number of events 


T T T T T T T T 

oooooooo 

cqcqcqcqcqcqcqcq 


u (— Dttjtt^ 

D~h+ ■ B^ ^ D-ii+Ui, 

D*-p+ 


yes 
yes 
yes 
yes 
no 
no 
no 
no 


(39.7 ±3.7) X 10-6 
(2.57 ± 0.10) X 10 
(2.94 ±0.19) X 10-^ 
(14.2 ± 1.4) X 10-3 


3377 ± 60 

1803 ± 43 
142 ± 13 
37 ± 16 
7.3 ±0.7 

43.8 ± 2.5 
281 ± 20 

34.5 ±4.6 


B+-. 


{neutral generic decays} 
■ {charged generic decays} 


no 
no 




114 ±7 
282 ± 11 



Table 2: Summary of B background modes taken into account for the likelihood model. When the 
yield is varied in the fit, the quoted number of events corresponds to the fit results. Otherwise the 
expected number, taking into account the branching ratios and efficiency, is given. 



are grouped into ten different classes that gather decays with similar kinematic and topological 
properties: nine for neutral B decays, one of which accounts for inclusive decays, and one for 
charged inclusive B decays. 

Table 2 summarizes the ten background classes that are used in the fit. When the yield of a 
class is varied in the Maximum Likelihood fit the quoted number of events corresponds to the fit 
results. For the other modes, the expected number of selected events is computed by multiplying 
the selection efficiency (estimated using MC simulated decays) by the branching fraction, scaled to 
the dataset luminosity (347 ft)-^). The world average branching ratios [3, 31] have been used for 
the experimentally known decay modes. 

4.2 THE MAXIMUM-LIKELIHOOD FIT 

We perform an unbinned extended maximum-likelihood fit to extract the inclusive B^ K^t:^t:~ 
event yield and the resonant amplitudes. The fit uses the variables mES, A£', the NN output and 
the SDP to discriminate signal from background. The At measurement allows the determination 
of mixing-induced CP violation and provides additional continuum background rejection. 

The selected on-resonance data sample is assumed to consist of signal, continuum background 
and B background components, separated by the flavor and tagging category of the tag side B 
decay. The signal likelihood consists of the sum of a correctly reconstructed ("truth-matched", 
TM) component and a misreconstructed ("self-cross-feed", SCF) component. 

The probability density function (PDF) 7-*? for an event i in tagging category c is the sum of 
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the probability densities of all components, namely 



7?? = AT . f'[ 



(1 ~ /sCF)'^sig-TM,i + /sCp'^sig-SCF,? 



^ggn 9tag,i^gg,tag) 'Pqq,i 



class 



Class 

+ ^B+jfB+3 2 ^tag,i^S+,tagj) ^S+.ij 



class 



+ E ^BOjrBojno,ij , (15) 

where A'^gig is the total number of K^tt^tt^ signal events in the data sample; f^^^ is the fraction of 
signal events that are tagged in category c; /scf is the fraction of SCF events in tagging category 
c, averaged over the DP; 'Pgig.xMi ^^'^ ^sig-SCFi products of PDFs of the discriminating 

variables used in tagging category c for TM and SCF events, respectively; N^^ is the number of 
continuum events that are tagged in category c; (j'tag.i is the tag flavor of the event, defined to be +1 
teg and -1 for a ^^^^g; A^^- tag 



for a and —1 for a ; A^^^tag parameterizes possible tag asymmetry in continuum events; 



"Pgg j is the continuum PDF for tagging category c; N^^^ (^dass) is the number of charged (neutral) 
-B-related background classes considered in the fit, namely one (nine); ^s+j {^B^j) is the number 
of expected events in the charged (neutral) B background class j; f%+j {f%Oj) is the fraction of 
charged (neutral) B background events of class j that are tagged in category c; A^+ ^^^ j describes 
a possible tag asymmetry in the charged B background class j; V^^ is the S+ background PDF 
for tagging category c and class j; and V^o -j is the neutral B background PDF for tagging category 
c and class j. Correlations between the tag and the position in the DP are absorbed in tag- flavor- 
dependent SDP PDFs that are used for charged B and continuum background. The PDFs Vx 
{X = {sig — TM, sig— SCF, qq, 5+, B^) are the product of the four PDFs of the discriminating 
variables ^, xi = ruES, X2 = AE, xs = NN output and the triplet X4 = {m', 0', At}: 

fe=l 

where i is the event index and j \s & B background class. The extended likelihood over all tagging 
categories is given by 

7 _ 

c=l i 

where A^'^ is the total number of events expected in category c. 

A total of 75 parameters are varied in the fit. They include the 12 inclusive yields (signal, four 
B background classes and seven continuum yields, one per tagging category) and 30 parameters 
for the complex amplitudes from Equation (11). Most of the parameters describing the continuum 
distributions are free in the fit. 



^Not all the PDFs depend on the tagging category. The general notations Px,i{j) and Vx.iij) used for simplicity. 
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4.2.1 THE At AND DALITZ PLOT PDFS 



The SDP PDFs require as input the DP dependent relative selection efficiency, e = e{m',6'), and 
SCF fraction, /scF = /scf("^', ^')- Both quantities are taken from MC simulation. Away from the 
DP corners the efficiency is uniform. It decreases when approaching the corners, where one of the 
three particles in the final state is close to rest so that the acceptance requirements on the particle 
reconstruction become restrictive. Combinatorial backgrounds and hence SCF fractions are large 
in the corners of the DP due to the presence of soft tracks. 
For an event i, we define the time-dependent SDP PDFs 

nig-TM,i(m', e', At) = £i{l- fscF,i) I det Ji\ \A^{At)\'' , (18) 

nig-scF,i(m', e', At) = Si fscF,i I det Ji\ \A^iAt)\'' , (19) 

where Psig-TM,i{fn' , 9' , At) and PsigscF, i{fn' , 9' , At) are normalized. The phase space integration 
involves the expectation values {e (1 — /scf) \det J\FkF*) and (e/scF |det J|FkF*) for TM and 
SCF events, where the indices k, a run over all resonances belonging to the signal model. The 
expectation values are model-dependent and are computed with the use of MC integration over the 
SDP: 

, fn e (1 - /scf) I det J| dmW , , 

e 1 - /scf) det J F,f; = ^° ,i 7 ' ' , ,1, - 20 

/o /o £ I det J| F^F* dm'dO' 

and similarly for (e/scF | det J\ Fk,F*), where all quantities in the integrands arc DP dependent. 

Equation (15) invokes the phase space-averaged SCF fraction /scf = (/sCF | det J| F^F*). The 
PDF normalization is decay-dynamics-dependent and is computed iteratively. We determine the 
average SCF fractions separately for each tagging category from MC simulation. 

The width of the dominant resonances is large compared to the mass resolution for TM events 
(about 8MeV/c^ core Gaussian resolution). We therefore neglect resolution effects in the TM 
model. Misreconstructed events have a poor mass resolution that strongly varies across the DP. It 
is described in the fit by a 2 x 2-dimensional resolution function 

RsC¥{m'r,9'^,mt,e't) , (21) 

which represents the probability to reconstruct at the coordinate {'m'j.,9'^) an event that has the 
true coordinate (m^,^^). It obeys the unitarity condition 

1 1 

J J i?sCF(m;, 9'„ m[, e[) dm',d9', = 1 , (22) 



and is convolved with the signal model. The i?sCF function is obtained from MC simulation. 

We use the signal model described in Section 2.1. It contains the dynamical information and is 
connected with At via the matrix element in Equation (11), which serves as the PDF. The PDF 
is diluted by the effects of mistagging and the limited vertex resolution [38]. The At resolution 
function for signal (both TM and SCF) and B background events is a sum of three Gaussian 
distributions, with parameters determined by a fit to fully reconstructed decays [41]. 

The SDP- and Ai-dependent PDFs factorize for the charged B background modes, but not 
necessarily for the neutral B background due to B^B^ mixing. 
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The charged B background contribution to the hkehhood (Equation (15)) involves the param- 
eter ^s+^tag! multiphed by the tag flavor gtag of the event. In the presence of significant tag- 
"charge" correlation (represented by an effective flavor-tag- versus-Dalitz-coordinate correlation), it 
parameterizes possible direct CP violation in these events. We use distinct SDP PDFs for each 
reconstructed B flavor tag, and a flavor-tag-averaged PDF for untagged events. The PDFs are 
obtained from MC simulation and are described by histograms. The At resolution parameters are 
determined by a flt to fully reconstructed 5+ decays. For the background class we adjust the 
effective lifetime to account for the misreconstruction of the event that modifles the nominal At 
resolution function. 

The neutral B background is parameterized with PDFs that depend on the flavor tag of the 
event. In the case of CP eigenstates, correlations between the flavor tag and the Dalitz coordinate 

are expected to be small. However, non-CP eigenstates, such as a^vr^, may exhibit such correla- 
tions. Both types of decays can have direct and mixing-induced CP violation. A third type of decay 
involves charged D mesons and does not exhibit mixing-induced CP violation, but usually has a 
strong correlation between the flavor tag and the DP coordinate (the D meson charge), because it 
consists of i?-flavor eigenstates. Direct CP violation is also possible in these decays, though it is set 
to zero in the nominal model. The DP PDFs are obtained from MC simulation and are described 
by histograms. For neutral B background, the signal At resolution model is assumed. 

The DP treatment of the continuum events is similar to that used for charged B background. 
The SDP PDF for continuum background is obtained from on-resonance events selected in the 
rriEs sidebands and corrected for feed-through from B decays. A large number of cross checks 
have been performed to ensure the high fidelity of the empirical shape used. Analytical models 
were found to be insufficient. The continuum At distribution is parameterized as the sum of three 
Gaussian distributions with common mean and three distinct widths that scale the At per-event 
error. This yields six shape parameters that are determined by the fit. The model is motivated 
by the observation that the At average is independent of its error, and that the At RMS depends 
linearly on the At error. 

4.2.2 DESCRIPTION OF THE OTHER VARIABLES 

The mES distribution of TM signal events is parameterized by a bifurcated Crystal Ball function [43- 
45], which is a combination of a one-sided Gaussian and a Crystal Ball function. The mean and two 
widths of this function are determined by the fit. A non-parametric function is used to describe the 
SCF signal component. The AE distribution of TM events is parameterized by a double Gaussian 
function. Misreconstructed events are described by a non-parametric function. 

Both TTT-ES and AE PDFs are described by non-parametric functions for all B background 
classes. Exceptions to this are the rriES PDFs for B^ D^tt^ and B^ J/ipK^ components, 
and the AE PDF for B^ £)-7r+, which are the same as the corresponding distributions of TM 
signal events. 

The rTT-ES and AE PDFs for continuum events are parameterized with an ARGUS shape func- 
tion [46] and a first-order polynomial, respectively, with parameters determined by the fit. 

We use non-parametric functions to empirically describe the distributions of the NN output 
found in the MC simulation for TM and SCF signal events, and for B background events. We 
distinguish tagging categories for TM signal events to account for differences observed in the shapes. 

The continuum NN distribution is parameterized by a third-order polynomial that is constrained 
to take positive values in the range populated by the data. The coefficients of the polynomial are 
determined by the fit. Continuum events exhibit a correlation between the DP coordinate and the 
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shape of the event that is exploited in the NN. To correct for residual effects, we introduce a linear 
dependence of the polynomial coefficients on the distance of the DP coordinate from the kinematic 
boundaries of the DP. The parameters describing this dependence are determined by the fit. 

5 FIT RESULTS 

The maximum-likelihood fit of 22525 candidates results in a — s- K^-k^-k' event yield of 2172 it 
70 and a continuum yield of 14272 ± 126, where the errors are statistical only. Figure 1 shows 

distributions of the likelihood ratio (signal/background) for all the events entering the fit and for 
the signal-like region. Figure 2 shows distributions of /S.E, mES) NN output, Ai/(T(Ai), where cj(Af) 
is the per-event error on Af, as well as the DP variables m' and 6' , which are enhanced in signal 
content by requirements on the signal-to-continuum likelihood ratios of the other discriminating 
variables. Figure 3 shows similar distributions for m^K^ii'^), m{K^7r^) and m{TT~^TT~). These 
distributions indicate the good quality of the fit. Signal enhanced distributions of At and At 
asymmetry for events in the regions of /o(980)i^g and (770) Kg are shown in Figure 4. 

In the fit, we measure directly the relative magnitudes and phases of the different components 
of the signal model. The results are given together with their statistical errors in Tabic 3. The 
measured relative amplitudes Ca, where the index represents an intermediate resonance, are used 
to extract the Q2B parameters, for which the definitions are given below. 

Table 3: Summary of fit results for the magnitudes \ca\ and phases (f) in degrees of the resonant 
amplitudes. The quoted error is statistical only. 



Resonance Name 


|Co-| 


(/> [degrees] 






(;^) [degrees] 


/o(980)i^^ 


4.0 


0.0 


2.8 ± 


0.7 


-88.6 ±21.3 


p0(770)K0 
/o(1300)KO 


0.10 ± 0.02 


58.6 ± 16.4 


0.09 ± 


0.02 


21.3 ±21.2 


1.9 ±0.4 


117.6 ±22.6 


1.1 ± 


0.3 


-15.2 ±23.8 


Nonresonant 


3.0 ±0.6 


13.8 ± 14.3 


3.7 ± 


0.5 


-16.2 ± 17.3 


is:*+(892)7r- 


0.136 ±0.021 


-60.7 ± 18.5 


0.113 ± 


0.018 


102.6 ±22.9 


ii:*+(1430)7r- 


4.9 ±0.7 


-82.4 ± 16.8 


7.1 ± 


0.9 


79.2 ± 20.5 


/2(1270)KO 


0.011 ± 0.004 


62.9 ± 23.3 


0.010 ± 


0.003 


-73.9 ± 27.8 




0.34 ±0.15 


68.7 ±31.1 


0.40 ± 


0.11 


154.5 ± 28.6 



For a resonant decay mode a which is a CP eigenstate, the following Q2B parameters are 
extracted: the angle 2/3efr, defined as: 

2/3eff(a) = argic^c*^) , (23) 
and the direct and mixing-induced CP asymmetries, defined as: 

^(-) = + ' (24) 
2Im{cacl) _ 

^^''^ = \r \2^\r |2 5 (25) 
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for a self-tagging resonant decay mode a such as iC*+(892)7r , the direct CP asymmetry is 

defined as: 



Ic-P - Ic l2 



where 



the relative isobar phase between two resonances, a and k is defined by 

A^((7, k) = arg{ccrcl) ; (27) 
and the similar quantity for a self-tagging resonant decay mode, a and its charge conjugate a is: 

A(l){a,a) = arg{c,^) ■ (28) 

recall that we use a convention in which the decay amplitudes have absorbed the phase from 
B'^B^ mixing. Finally, we also extract the relative fraction / of a Q2B channel a, which is calculated 
as: 

^ {\c,\^ + \c,\^){F,F*) 

{F^F:) = j F^F:ds+ds- . (30) 

Two approaches were used to extract the statistical uncertainties of Q2B parameters. The first 
approach uses a linear approximation, whereby errors are evaluated assuming that the likelihood 
function is a multivariate Gaussian, which is defined by the fit result and covariance matrix. In 
the second approach we perform a likelihood scan, fixing the scanned Q2B parameter at several 
consecutive values and repeating the fit to the data. The error on the Q2B parameter is obtained 
from the shape of the likelihood function near the minimum. Since the Q2B parameters are not 
directly used in the fit, we fix instead certain parameters that allow the resulting likelihood curve 
to be trivially interpreted in terms of the Q2B parameter of interest. The likelihood scan approach 
does not rely on any assumption about the covariance matrix. We use this approach wherever the 
fit variables that determine a Q2B parameter are found to have non-Gaussian errors. 

The Q2B parameters and fit fractions are given together with their statistical and systematic 
errors in Table 4. It is indicated in the table whether the statistical error has been computed 
by a likelihood scan or using the linear approximation. Systematic uncertainties are discussed in 
Section 6. Results of likelihood scans are shown in Figure 5 in terms of = — 21n(>C). 

The measured values of 2/3cff(/o(980)E:^) and 2/3cfr(/5°(770)Kg) are both consistent with the 
SM predictions. For the former, the measured value is higher by 2.1 standard deviations compared 
to that for b ccs. This is unlike the tendency of other results in 6 — > qqs transitions. In 
addition to this, 2/3es{fo{^80)Kg) = is excluded at 4.3(7 significance. The phase difference 
A(j){K*~^ {892)tt~ , K*~ {892)Tr~^) is measured here for the first time. A mirror solution at ~ 35° is 
disfavored at 3.7o- significance. The interval -102.0° < A(j}{K*+ {892)^- , K*' {892)^+) < 135.7° is 
excluded at 95% confidence level (CL). 

Note that the values of the Q2B parameters S and C that are obtained from our fit variables 
must take values within the physical boundary S"^ + < 1, in contrast to the values obtained 
from Q2B analysis. The value of S{f o{980) Kg) that is obtained is close to the physical boundary 
and consequently has a highly non-Gaussian uncertainty, as shown in terms of in Figure 5. The 
positive uncertainties on S{fo{980)Kg), including systematics, extracted from this distribution 
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Table 4: Summary of results for the Q2B parameters. The first quoted error is statistical, the second 
is systematic and the third is DP signal model uncertainty. Parameters for which the statistical 
error have been obtained from, a likelihood scan are mxirked by f. Phases are in degrees and relative 



fractions i/n %. 



Parameter 


Value 


Parameter 


Value 


t2/3eff(/o(980)Kg) 
tS(/o(980)/i'0) 

/(/o(980)i^0) 


U.oO±U.Z( ztU.U / ±U.U4 

(89+20 ± 5 ± 8)° 

-n 04+0.07+0.05 _L Q no 
^•^^-0. 02-0.03 ^-^^ 

14.3iij±1.5±0.6 


C(p (nU)Kg) 
t2/?eff(p°(770)ii:g) 

/(/(770)ifg) 


n no -L n ot -l n no -l n 

(37li^ ± 5 ± 6)° 
0.6ltg:^^± 0.09 ±0.08 
9.0±'l.4±l.l±l.l 


Acp{K*+{m2)TT-) 
tA(/)(if*(892)7r) " 
/(i^*(892)^) 


-0.18 ±0.10 ±0.03 ±0.03 
(-164 ±24 ±12 ±15)° 
11.7± 1.3± 1.3±0.6 


tA<^(/oi^°,pOi^°) 


(-59tJ?±6±6)° 


/(X*(1430)7r) 
/(/o(1300)ifg) 
f{Xco{lP)K%) 


38.9±2.5±0.7±1.3 
6.3 ± 1.3 ±0.6 ±0.3 
1.2 ±0.5 ±0.0 ±0.1 


./(iVi?) 
/(/2(1270)i^g) 


25.6 ±2.5 ±1.9 ±0.5 
2.1 ±0.8 ±0.0 ±0.2 



"Abbreviation for A4>{K*+{892)tt- , K'- {892)tt+). 



are 0.09 at the 32% CL and 0.31 at the 5% CL. The systematic error on 5" is derived from the 
corresponding errors on 2/3efr and C. 

To vahdatc the presence of the fx{1300)Kg, f2{1270)Kg and Xco{^P)Kg resonant modes in 
the signal, we performed fits with three reduced signal models, in which we removed these modes 
one by one. The differences between these fits and the nominal one in terms of for the models 
missing /x(1300)K§, /2(1270)K§ and Xco{'i-P)K^ are 74, 37 and 35, respectively. 

As a validation of our treatment of the time-dependence we allow r^o and Amd to vary in 
the fit. We find r^o = (1.579 ± 0.061) ps and Anid = (0.497 ± 0.035) ps"^ while the remaining 
free parameters are consistent with the nominal fit. In addition, we performed a fit floating the S 
parameters for S° J/i^K^ and S° ijj{2S)K^ events. We measure S = sin(2/?) = 0.690 ±0.077 
and 0.73 ± 0.27 for J/tpK'^ and ip{2S)K'^ respectively. These numbers are in agreement with 
the current world average for sin(2/3). To validate the SCF modeling, we leave the average SCF 
fractions per tagging category free to vary in the fit and find results that are consistent with the 
MC estimation. 



6 SYSTEMATIC STUDIES 

The contributions to the systematic error on the signal parameters are summarized in Table 5. 

To estimate the contribution to K^tt+tt" decay via other resonances, we have first fitted 

the data including these other decays in the fit model. We considered possible resonances, including 
tj(782), p°(1450), /9°(1700), /o(1710), /2(1810), K*^{1680), ^2=^(1430), Xc2{lP) and a low mass 
S wave a. A RBW lineshape has been used to parameterize these additional resonances, with 
masses and widths from [31]. As a second step we have simulated high statistic samples of events, 
using a model based on the previous fits, including the additional resonances. Finally, we fitted 
these simulated samples using the nominal signal model. The systematic effect (contained in the 
"DP model" field in Table 5) is estimated by observing the difference between the generated values 
and the fitted values corresponding to the generated samples. This systematic effect is quoted 
separately. 
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Figure 1: Distributions of the likelihood ratio (signal/background) for all events entering the fit 
(left) and in the signal-like region (right). The dots with error bars give the on-resonance data. 
The solid histogram shows the projection of the fit result. The dark, medium and light shaded areas 
represent respectively the contribution from continuum events, the sum of continuum events and the 
B background expectation, and the sum of these and the misreconstructed signal events. The last 
contribution is hardly visible due to its small fraction. In both distributions the D~ir^ and J/tpKg 
bands are removed from the DP. 

We vary the mass, width and other parameters (if any) of all components in the fit within their 
errors, as quoted in Table 1, and assign the observed differences in the measured amplitudes as 
systematic uncertainties ("lineshape parameters" in Table 5). 

To validate the fitting tool, we perform fits on large MC samples with the measured proportions 
of signal, continuum and B background events. No significant biases are observed in these fits. The 
statistical uncertainties on the fit parameters are taken as systematic uncertainties ("Fit bias" in 
Table 5). 

Another major source of systematic uncertainty is the B background model. The expected 
event yields from the background modes are varied according to the uncertainties in the measured 
or estimated branching fractions. Since B background modes may exhibit CP violation, the cor- 
responding parameters are varied within their uncertainties, or if unknown, within the physical 
range. As is done for the signal PDFs, we vary the At resolution parameters and the flavor-tagging 
parameters within their uncertainties and assign the differences observed in the data fit with respect 
to the nominal fit as systematic errors. The systematic uncertainties from these sources are listed 
as Background " in Table 5. 

Other systematic effects are much less important for the measurements of the amplitudes and 
are combined in the "Other" field in Table 5. Details are given below. 

The parameters of continuum PDFs are determined by the fit. No additional systematic un- 
certainties are assigned to them. An exception to this is the DP PDF: to estimate the systematic 
uncertainty from the m-^s sideband extrapolation, we select large samples of off-resonance data by 
loosening the requirements on AE and the NN output. We compare the distributions of m' and 0' 
between the m-^s sideband and the signal region. No significant differences are found. We assign 
as systematic error the effect seen when weighting the continuum DP PDF by the ratio of both 
data sets. This effect is mostly statistical in origin. 

The uncertainties associated with Arrid and r are estimated by varying these parameters within 
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Figure 2: Distributions of (top to bottom, left to right) AE, m-Es, NN output, At/a{At), m' and 
9' for samples enhanced in — > K^tt'^tt^ signal. The dots with error bars give the on-resonance 
data. The solid histogram shows the projection of the fit result. The dark, medium and light shaded 
areas represent respectively the contribution from continuum events, the sum of continuum events 
and the B background expectation, and the sum of these and the misreconstructed signal events. 
The last contribution is hardly visible due to its small fraction. In all these distributions the D^n^ 
and J/ip bands are removed from the DP. 
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Figure 3: Distributions of m{K^7r~^) (top left), m{K^TT~) (top right), and m['K~^7r~) (bottom) 
for samples enhanced in — > K^t^^tt^ signal. The solid histogram shows the projection of the 
fit result. The dark, medium and light shaded areas represent respectively the contribution from 
continuum events, the sum of continuum events and the B background expectation, and the sum of 
these and the misreconstructed signal events. The last contribution is hardly visible due to its small 
fraction. In all these distributions the 0"-^^ and J/ipK'^ bands are removed from the DP. 
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Figure 4: Distributions of At when the i?tag ^-^ ^ -^'^ (^^p), (middle) and the derived At asymme- 
try (bottom). Plots on the left (right) hand side, correspond to events in the /o(980)-fr^ (p^{nQ)K^) 
region. The solid line is the total PDF, the dashed line is the continuum only PDF and points with 
error bars represent data. These distributions correspond to samples enhanced in K^'k^-k~ 
signal, where the D~Tr^ and .J/ip Kg bands are removed from the DP. 
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Figure 5: Results of likelihood scans in terms of = — 21n(£) for the Q2B parameters 
(left to right, top to bottom): S{fo{980)K°), S{p°{770)K^), 2/3eff(/o(980)Kg), 2/3eff(p°(770)Kg), 
A0(/°(98O)Kg,p(77O)Kg) and A(/)(K*+(892)7r-, i^*-(892)7r+). The solid (dotted) curves corre- 
sponds to the total (statistical) error. Indicated by solid (dashed) rectangles are the parameter 
values corresponding to la (2a). The arrows mark the measured values in b ^ ccs transitions [3]. 
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the uncertainties on the world average [31]. 

The signal PDFs for the At resolution and tagging fractions are determined from fits to a control 
sample of fully reconstructed B decays to exclusive final states with charm, and the uncertainties 
are obtained by varying the parameters within the statistical uncertainties. 

The average fraction of misreconstructed signal events predicted by the MC simulation has been 
verified with fully reconstructed B Dp events [38]. No significant differences between data and 
the simulation were found. We vary /scf for all tagging categories relatively by 25% to estimate 
the systematic uncertainty. Tagging efficiencies, dilutions and biases for signal events are varied 
within their experimental uncertainties. 

7 SUMMARY 

We have presented preliminary results from a time-dependent Dalitz plot analysis of B^ Kgir^iT^ 
decays obtained from a data sample of 383 million T{AS) BB decays. We measure 15 pairs of 
relative phases and magnitudes for the different resonances, taking advantage of the interference be- 
tween them in the Dalitz plot. Prom the measured decay amplitudes, we derive the Q2B parameters 
of the resonant decay modes. In particular, the mixing- induced CP asymmetry S is extracted from 
the measured amplitudes. The measured values of 2/3cff in B^ decays to fo{980)Kg and (770) Kg 
are (89^20 =t 5 it 8)° and (37^^^ ± 5 ± 6)°, respectively. These results are both consistent with the 
SM predictions, but in the case of B^ fo{980)Kg the measured value is higher by 2.1 standard 
deviations compared to that for b ccs. This is unlike the tendency of other results in 6 — ^ qqs 
transitions. Also, 2/3cfT(/o(980)X^) = is excluded at 4.3 o" signficance. 

In decays to A'*(892)7r we find Acp = -0.18 ± 0.10 ± 0.03 ± 0.03. The phase difference Ac/) 
between the amplitudes of S° — i^*+(892)7r~ and B^ K*~(892)7r+ is measured for the first time. 
We find = (—164 it 24 it 12 it 15)°, and mirror solutions are disfavored at ^ 3.7a significance. 
The interval —102° < Acp < 136° is excluded at 95% confidence level. Our results may be used to 
extract the CKM angle 7 following the methods proposed in Refs. [18-21]. 
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C{fo(.980)Kl) 


2/3eff(/o(980)ifO) 


/(/o(980)ii:0)[%] 


DP model 


0.04 


7.8 


0.6 


Lineshape parameters 


0.06 


3.9 


1.0 


B background 


0.03 


2.7 


0.4 


Fit bias 


0.01 


1.1 


1.0 


Other 


0.03 


1.4 


0.1 


Sum w/o DP model 


0.07 


5.1 


1.5 


Total sum 


0.08 


9.3 


1.6 






C{p^{770)K's) 


2/3eff(p°(770)i^O) 


/(pO(770)i^O)[%] 


DP model 


0.06 


5.9 


1.1 


Lineshape parameters 


0.04 


3.6 


0.3 


B background 


0.06 


3.8 


0.1 


Fit bias 


0.02 


0.4 


1.0 


Other 


0.02 


1.0 


0.1 


Slim \\-/o DP ukkU'I 


0.08 


5.3 


1.1 


Total sum 


0.10 


7.9 


1.5 





Acp{K*+7r-) 




[%] A(l){K*+7r-,K*-Tr+) 




DP model 


0.03 


0.6 


15.0 


6.0 


Lineshape parameters 


0.01 


0.2 


4.3 


4.2 


B background 


0.03 


0.3 


4.5 


4.3 


Fit bias 


0.01 


1.2 


9.7 


0.3 


Other 


0.00 


0.1 


2.6 


1.7 


Sum w/o DP model 


0.03 


1.3 


11.8 


6.3 


Total sum 


0.05 


1.4 


19.1 


8.7 



Table 5: Summary of systematic uncertainties. Errors on 2/?eff and A(f) are given in degrees and 
relative fractions in %. K*^ refer to K*^{892). 
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